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Abstract 

The temperature dependence of the elastic scattering processes qq' 
qq' and qq' qq', with q, q' = u,d,s is studied as a function of the 
scattering angle and the center of mass energy of the collision within 
the framework of the SUf{3) Nambu-Jona-Lasinio model. Critical 
scattering at threshold is observed in the qq' — > qq' process, leading to 
an enhancement of the cross section as occurs in the phenomenon of 
critical opalescence. Transport properties such as viscosity, mean free 
paths and thermal relaxation times are calculated. Strangeness en- 
hancement is investigated via the chemical relaxation times, which are 
found to be considerably higher than those calculated via perturbative 
QCD. A comparison with the experimental values for the strangeness 
enhancement in 5 + 5" collisions leads to an upper limit of 4 fm/c for 
the lifetime of the plasma. 
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I. INTRODUCTION 



With the proposed experiments for heavy ion coUisions at the LHC coming closer, 
the onus has settled on the theoretical physics community to provide a clear, coherent 
and credible description of such collisions. One of the major goals, that of detecting 
the presence of plasma formation of the constituent quarks and gluons is for both 
theorists and experimentalists alike, a formidable task. Such a phase is believed to 
be formed as matter is deconfined. 

From the theoretical side, a description that goes beyond the standard equilib- 
rium statistical approaches is necessitated by the short collision time of the ions, 
since this may preclude equilibration of all degrees of freedom. Thus a transport 
theory that is consistently based on quantum chromo dynamics (QCD) itself is re- 
quired. While such formulations have been set up by some authors |jl|J^, they have 
not reached the stage where they can be implemented to provide original physical 
results. 

One step in this direction is to examine the consequences of a transport theory 
based on an effective model of QCD, that exhibits some important features of the 
theory. In particular, it is desirable to see if any major effects due to the exhibited 
feature are observable in the final analysis. 

It is well-known that the low energy particle sector is well described by effective 
chiral theories of QCD. We therefore choose to implement the Nambu-Jona-Lasinio 
(NJL) model in its three flavor extension for such a study. This model has the 
advantage that, being a model, it can be studied over the entire temperature range. 
The NJL model also offers a simple intuitive view of chiral symmetry breakdown and 
restoration via a realization of the pairing mechanism of quark and antiquark states 
of an effective interaction in a similar sense to the BCS theory of superconductivity. 
The model suffers however from several unaesthetic diseases: the simple form of a 
point-like interaction, while making the model analytically tractable, also ensures 
its nonrenomalizability, and a cutoff scale A for the theory must be introduced. A 
confinement / deconfinement phase transition cannot occur. Rather a "soft" form of 
deconfinement is observed at the so-called Mott temperature, at which the mesons 
dissociate into their constituents 0. 

The NJL model, being a strong coupling field theory, cannot be handled in 
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perturbation theory. The by now standard technique is a self-consistent expansion 
of quantities to be calculated in the inverse number of colors l/A'^c 0. It is this 
expansion, taken to lowest order, that we use in our model. 

In order to make some progress in the ambitious task presented, we have already 
studied hadronization processes of quarks and antiquarks into two mesons, qq' 
MM' P], within the SUf{3) NJL model. In addition, we require the elastic scattering 
cross sections for both quarks and quarks qq' — * qq' and quarks and antiquarks qq' —>■ 
qq' within the three flavor model. The cross sections for the third possible class, 
qq' qq', follow immediately from qq' —>■ qq' by time reversal invariance, as long as 
no chemical potential is involved. It is the purpose of this paper to calculate these 
cross sections, together with the attendant transport properties (mean free paths and 
shear viscosities) in the SUf{3) model for realistic values of the current quark masses 
rriQu = fnod = 5.5 MeV, = 140 MeV. We note that some elastic scattering cross 
sections have been calculated in the two flavor sector with the additional restriction 
of the chiral limit condition = mod = in Ref. P]. As will be seen, the extension 
presented here indicates significant increases in the cross sections due to additional 
available channels. 

The quark quark scattering cross sections are relatively featureless. On the 
other hand, quark antiquark scattering shows a threshold divergence at the Mott 
temperature Tmtt, at which the pion dissociates into its constituents. At higher 
temperatures that are to be expected in the plasma phase, a large resonant structure 
is observed, that is a remnant of this feature of critical scattering. This is however 
not directly observable. When these cross sections are included in a fully consistent 
transport theory based on this chiral model of QCD, it is hoped that this may lead 
to observable effects. 

We have also studied transport coefficients witin this model. We find that the 
viscosity coefficient rj/T^ lies in the range 0.83 to 1.0 per flavor degree of freedom 
for light quarks, and is slightly lower for strange quarks, being 0.68 to 0.75, i.e. 
the total viscosity is found to lie between 3.2 and 4.0 considering only two flavors, 
while for three flavors the total viscosity would lie between 4.6 and 5.5. Perturbative 
QCD estimates lie slightly higher than this value: a variational calculation ||10[ (see 
also ||TT]]) gives rj/T^ ^ — 1.16/(a;^logas) for two quark flavors, where as is the 
strong coupling strength. One flnds the value rj/T^ ^ 6.3 at = 0.6. This should 
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be viewed somewhat critically, since this perturbative formula is strictly valid for 
Q^s <^ 1, and is an extremely sensitive function of Ug- In the strongly interacting 
non-perturbative regime, it is not applicable. Hydrodynamical estimates, on the 
other hand |12[, give a wide range for this quantity, 2 < rj/T^ < 3tT, where r is the 
proper time for the expansion of the system. We also calculate mean free paths for 
light and heavy flavors and these are found to be of the order of magnitude of 1 fm. 
The thermal relaxation times are found to vary between 1-2 fm/c, with that of the 
strange quark being slightly higher. 

In order to discuss strangeness enhancement, we examine the chemical relaxation 
time for strange quarks within the model. These, lying between 11 fm/c and 17 fm/c, 
are considerably higher than the naive perturbative QCD prediction of Ref. [0 that 
proceeds via gluon fusion. The origin of this discrepancy can be traced back to the 



value of the strange quark mass m^. In the perturbative QCD calculation of Ref. |]13 



the physical mass was assumed to be equal to the current quark mass value. In the 
NJL model, on the other hand, a significantly enhanced value of the strange quark 
mass is still found at temperatures of the order of 250 MeV ||^. In fact it is at 
least twice the current quark mass value. A recalculation of the perturbative QCD 
calculation using the temperature dependent NJL masses gives relaxation times of 
the same order of magnitude as our calculation. We note that an assessment of 
the relevance of quark antiquark scattering for the strangeness enhancement that is 
observed e.g. in S* + collisions, depends on the plasma lifetime. Our calculations 
would attribute the observed enhancement to these processes for a plasma lifetime of 
ca. 4 fm/c. Such a value is supported by hydrodynamical calculations [^^-|TB[. This 
places an upper limit on the plasma lifetime. For shorter times, other mechanisms, 
including hadronization and final state interactions must also be invoked. 

Our paper is organized as follows: In Section 0, we first give a brief review 
of essential formalism, then we discuss quark quark scattering in the nonstrange 
and strange sectors. Following this, quark antiquark scattering for both nonstrange 
and strange partons is presented. Transport properties are derived and discussed 



in Section |ITT[ The possibility of strangeness enhancement due to quark antiquark 
scattering is discussed in Section 0. A general discussion of our methods and results 
is presented in Section |V]. We summarize and conclude in Section [VJ. 
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II. SCATTERING PROCESSES 



A. Review of Formalism 

The model which we employ for our calculations is the SUf{3) version of the 
Nambu-Jona-Lasinio (NJL) model considered at finite temperature. In this 
subsection, we give a brief review of the basic formalism and define our notation, 
without giving derivations of our formulae. For a more detailed treatment, the 
reader is referred to Refs. P,p|,p|,p^ . The Lagrangian for this model is 

C= ^/(^^-mo/)^/ + Gf:[(^A»2 + (^^75A»2] (2.1) 

f=u,d,s a=0 

-/s:[det^/^(l + 75)?/' + det^/'(l-75)^/^] , 

where G and K are coupling constants with dimensions [MeV]^^ and [MeV]~^ respec- 
tively, and niQf are the current masses for quarks of flavor /, which explicitly break 
chiral symmetry. The determinantal term leads to the appearance of six fermion 
vertices, which, in a mean field approximation, can be reduced to an effective four 
point interaction by contracting out ipip pairs Under the additional assumption 
of SUf{2) isospin symmetry, m„ = m^, this leads to the effective Lagrangian 

8 



J2 

f=u,d,s 



a=0 



(2.2) 



+ K,, [(^AV)(V^AV) + (v^aV)(^aV) 
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(^«75AV)(^n5AV) + (V^^75AV)(V'n5AV) 



with the effective coupling constants 



Kf = G± \Kg' 



Kt = Kf = G± iK^" 



J\l — ^^-2 — --3 — ^ ^ 2 

G ± |K(4^" - 



±^K{g'' - g^ 



(2.3) 



expressed in terms of the couplings G, K, and the trace of the Green function 
Sf{x, x) via 



gf = NjtT^Sfi 



X, X 



47r2 



mfA{mf,iif) 



(2.4) 
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where, in the case of finite temperature and chemical potential, 

A(m/,/ij) = ^ Ve*""^ / J— 7^ 2T • (2-5) 

|p1<a 

In this expression, ujn are the Matsubara frequencies for fermions, ujn = (2n + l)7r//5, 
and the sum runs over all positive and negative integer values. Note that this 
quantity is derived from the finite temperature Green function 

' (3^ J (27r)3 7o(za;„ + /i;)-7p-mj ^ ' 

Explicit expressions for A(m, /x), which are suitable for numerical evaluation, can 



be found in Ref. |T7]. Note that we employ a three momentum cutoff A in order 
to make the integral in Eq. ( |2.5| ) finite. The physical quark masses m^, are 
calculated from the three coupled gap equations 



rrii = niQi —miA{mi, /ij) + mjA{mj,iXj)mkA{mk, fXk) , (2.7) 



where i, j, k are three pairwise distinct fiavors. An explicit numerical solution of 
Eq. can be found in Ref. [|[. 

Mesons are constructed by calculating the quark-antiquark effective interaction 
within the random phase approximation 0,^. In the pseudoscalar sector, this leads 
to the explicit form for the pion and kaon propagators 



2K^ 

The irreducible pseudoscalar polarization function II^j, required in Eqs. (|2.8| ) and 
( p.9|) is given explicitly by |^ 



n;/,(po,p) = -^{Aimj,fif) + Aimj,,fif,) (2.10) 
where Bq is the analytical continuation of 
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Bq{p, rrif, fif, mf',fif>,iu^) = (2.11) 



13 ^ J (27r)3 [(z^„ + fi^y - Ej] [{tu^ - ivra + ^lf'f - E}] 



{Ef = ^(p + mj, Efi = \J{(i — p)"^ + rn^fi) to the real axis. Explicit expressions for 
Bq can again be found in [|17|. 

For the r] and the situation is somewhat more involved due to the mixing terms 
in Eq. ( p.2|) . The quark-antiquark scattering matrix containing these particles is a 
2x2 matrix 



with 



AC-B' 



A = K^-^ICinuu + ^^s,) (2.13a) 



B = i^o+ + ^muu - nf,) (2.13b) 

c = ir+ - ^/c(2n^, + nf,) (2.13c) 

IC = K^Ki-K',, . (2.13d) 

Since we do not require the meson masses and their static properties explicitly in 
this work, we do not display the technical details needed for this here, but rather 
refer our reader to Ref. 

In the scalar sector, the NJL model in its SUf{3) version contains nine mesons 
(cr^o, an±, (Tko, ctj^, aK±, cr, a'), that accompany the nine mesons (vr^, Tr^*", K^, 

, 7], 7]') in the pseudoscalar sector. The propagators for the scalar particles can be 
immediately obtained from Eqs. ( |2.8| )-( prT3D by replacing the pseudoscalar coupling 
strengths by the scalar ones K~ and replacing the pseudoscalar polarization by 
the scalar polarization 

nf/,(po,p) = -^[A{mf,iif) + A{mf,,iif.) (2.14) 

+ \{mf + mfif - {po + l^f - jJ^^y Bo{p,mf,iJ,f,mf,iJf,po)} . 

Via Eqs. (|2.8|), ( |2.9| ) and ( |2.12| ), the masses of the pseudoscalar mesons as well as 
those of the scalar mesons with the appropriate changes mentioned above, have been 
obtained as a function of temperature. For explicit values see Ref. M. 
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B. Quark Quark Scattering 

In this subsection, we classify and illustrate via calculation the possible inde- 
pendent elastic scattering collision processes within the SUf{3) flavor quark-quark 
combinations. 

For elastic scattering processes of the type qq' qq' , via simple combinatorics 
we have six different possibilities available. However, due to isospin symmetry, the 
processes us us and ds ds lead to the same cross section. The same holds 
true for the processes uu uu and dd — dd, so that we have in total only four 
independent processes for quark-quark scattering. These are listed in Table |. The 
Feynman diagrams appropriate for such scattering in lowest order in terms of a 
1/iVc expansion, where is the number of colors, are shown in Fig. |l|. Here both 
t-channel and w-channel diagrams can occur, with the species of the exchanged 
meson depending on the specific process in question. These are also given explicitly 
in Table | for each of the independent processes. We note that there is at least one 
scalar and one pseudoscalar exchange channel for each diagram. In this point, this 
SUf{3) chiral model differs markedly from the SUf{2) chiral model, as is evident 
from the fact that the SUf{2) model supports only one scalar meson. This will be 
seen explicitly in the numerical calculations. 

In the following, we will derive as far as possible general expressions for the 
cross sections. We illustrate these with explicit calculations for the specific process 
uu —>■ uu. 



1. Analytical Calculations 

The matrix elements corresponding to the Feynman graphs of Fig. ^ can in 
general be written as 



iMt = Sc-,,c3Sc2,cMP3)Tu{pi) iVfi^pi - ps) u{pi)Tu{p2) 

+ ^ci,c3^c2,cMP3)i^l5T)u{pi) iV^{pi-p3) u{pi){i-i^T)u{p2) (2.15) 



iMu = 6c-,,c45c2,c3uiP4)Tu{pi) liV^i^pi - p^)] u{ps)Tu{p2) 

+ 6c,,c^6c,,oMP4)iil5T)u{pi) \iV^{pi-p4)] u{p^){i-i^T)u{p2) 



(2.16) 
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Here TVfT is a symbolic expression for the sum over all scalar t exchange channels of 
flavor factors times particle propagator, with TV[T being the same for pseudoscalar 
t graphs. The symbols TV^T and TV^T are the analogous quantities for the u- 
channel. Note that the matrix elements ( |2.15| ) and ( |2.16| ) describe the scattering 
due to the exchange of colorless mesons. The scattering cross sections due to the 
exchange of color octet states can be estimated to differ from Eqs. ( p.l5|) , ( p.l6| ) by 
a factor 4/9 and are neglected here. 

To calculate the cross sections, we require the squared matrix elements, summed 
over spin and color degrees of freedom and averaged over the incoming states, 

^ Y,\Mt-Muf . (2.17) 



4Ar2 

c s,c 

General expressions for these, apart from the flavor factors stemming from the Gell- 
Mann matrices, are given in Appendix ^. The flavor factors, together with the 
propagators of the exchanged mesons, have to be specifled for each individual pro- 
cess. We illustrate this procedure explicitly here for the process uu uu. From 
Table |, we read off that the t-channel exchange proceeds via an 7r°, rj and rj' in the 
pseudoscalar sector and via a (77^0, a and a' in the scalar sector. This means that 
the propagators for the t-channel have the form 



f2.18) 



where Vt, is deflned in Eq. ( p.8|) and A, B, C and /C are given by Eq. ( p.l3|) . The 
propagator for the scalar t-channel has the same form with, however, the coupling 
constants and polarization functions suitably altered, as was discussed in Sec. |11 A| . 
The M-channel exchange proceeds via the exchange of the same particles. Thus for 
Vu we have the same form as in Eq. ( p. 181) . 

In calculating the cross section, we conflne ourselves to the situation in which 
the center of mass system of the incoming particles is at rest relative to the medium. 
In addition, in this paper, we will examine the cross sections at flnite temperature, 
and in what follows, set the chemical potential to zero. In this case the differential 
cross section can be written as 

da 1 1 ■s-^ , . , . , ,2 



dt 167rs(s - 4m2) 4Ar2 



Y^\Mt-MuV ■ (2.19) 



For the total cross section, we include a Fermi blocking factor for the flnal states 
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ldt^[l-MPV~s/2)]' , (2.20) 



cr 



with the Fermi function fpix) = l/(expx + 1). Note that the total cross section 
depends only on the temperature T and the center of mass energy ^/s 0. 

The process uu — uu is an example of the elastic scattering of quarks that is 
already present in the SUf{2) model. It differs from this three flavor calculation in 
that the additional exchange channels 0",^, rj, rj' and a' are missing in that case. For 
a pure SUf{2) calculation, only the terms 

TVfT = =- (2.21a) 

i-4Gns,(o,v^) ^ ^ 

TVfT = — =- (2.21b) 

TV^T = .^.^p^r, (2.21c) 
1 - 4Gn^^(0, v-m) 

occur, leading us to the result of Ref. 0. The appearance of the new exchange 
channels in both pseudoscalar and scalar sectors leads to a significant enhancement 
in the numerical values of the cross sections, as will be shown below. 

One notes from Table | that the ud — ud cross section follows via exchange of 
the same scalar and pseudoscalar mesons in the t-channel, but a restricted set in the 
^-channel is admitted. The elastic ss ss processes couple again to the same set of 
scalar and pseudoscalar mesons in both t- and w-channels. In this case, no coupling 
to the TT or o"^ can occur. When one light and one heavy quark are scattered, one 
finds that virtual K and o"x exchanges are possible in the u-channel. This leads to 
some technical difficulties that are due to the kinematics. One would expect that the 
kaon (and a^) propagator should not depend on the sign of the zero component of 
the exchanged momentum. This corresponds, via Eqs. ( p.8|) , (p.9|) , to the condition 
ngs(po,p) = ^qs{—po,P), or, through Eqs. (|2.10|) , ( p.l4| ), to the equivalent statement 



Bo{p,mq,fiq,rns,fis,Po)=Bo{p,mg,fiq,ms,fis,-po) ■ (2.22) 

However, this symmetry is lost in a straightforward evaluation which is an 
artifact of our regularization procedure. This can be explicitly traced back to the 
fact that a shift in variable is performed in certain terms occurring in Ilgj. Such 
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variable shifts are strictly non admissable for divergent integrals, but since they 
contribute only to order 1/A, they are usually ignored. The symmetry condition 
( p.22| ) is in any event still fulfilled if the constituent quark masses are equal (i.e. for 
pions, etas and their scalar partners) or the three momentum argument vanishes. For 
evaluations in the kaonic sector, however, the violation of ( p^.22[ ) becomes evident. 
Accompanying this artifact is the further cutoff and regularization scheme induced 
artifact that the imaginary part of the polarization is discontinous at po = 0, i.e. at 
the kinematical point that is required e.g. for the process uu ss in the center of 
mass system. We resolve this problem by hand by replacing -Bo(p, ^q, fJ'q, ^s, fJ^s,Po) 
with the symmetric form 

^ [Bo{p,mg,^g,m„fi„po) + Bo{p,mg, fig,m„ ^„ -po)] (2.23) 

whenever dealing with kaons. 



2. Numerical Results 



Our numerical calculations were performed using the parameter set niQu = rriQd = 
5.5 MeV, mo, = 140.7 MeV, GA^ = 1.835, A^A^ = 12.36 and A = 602.3 MeV. This is 
the same parameter set that was used in Ref . [§ , so the numerical results concerning 
the static mesonic properties can be obtained from this reference. In particular it 
was demonstrated in Ref. [^] that at the pionic Mott temperature Tmtt = 212 MeV, 
the pion mass becomes equal to the mass of its constituents m^r = 2m„ and the pion 
becomes a resonant state. The same happens with the kaon at T]\jk = 210 MeV, 
the r] at Tmtj = 180 MeV, and the a at Tmo- = 165 MeV. At these temperatures, the 
respective particles become unbound. This effect in a rather crude fashion models 
the deconfinement transition within the NJL model 0. 

We first compare our results for three flavors with the corresponding calculation 
of Ref. 0. In Fig. ^ we show the total cross section for the process uu uu at 
T = 215 MeV as a function of ^/s. The temperature chosen lies slightly higher than 
the pion Mott temperature Tmtt] we are therefore in the plasma phase, in which the 
model may be regarded as physically realistic. As in Ref. [§, our calculations are 
shown for center of mass energies < 2^Ja^~+^v^, i.e. they are restricted by the 
natural cutoff of this model. Numerically we find that the SUf{3) calculation yields 
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a far larger cross section for uu uu than the corresponding two flavor case, the 
difference being a factor of 3-4. Thus the greater number of exchange channels that 
become available for three flavors of quarks significantly enhances these scattering 
processes, and the SUf{2) result is only recovered on directly eliminating these. The 
remainder of our calculations are shown for the three flavor sector only. 

In Fig. ^, we show the total qq scattering cross sections that involve nonstrange 
quarks as a function of ^/s and for two values of the temperature, T = 215 MeV 
and T = 250 MeV. The curves for both uu —>■ uu and ud — > ud are essentially 
flat and display no particular structures. The ud ud scattering cross section lies 
at a given temperature slightly higher than the uu uu cross section, with both 
between 1.2 and 1.6 mb for T = 215 MeV. This difference can be attributed to 
the different flavor factors accompanying the various mesonic states. Note that the 
ud scattering in fact has less exchange mesons available in the w-channel, since no 
neutral particles are admissible in this channel. As the temperature is increased 
from T = 215 MeV to T = 250 MeV, the cross sections become somewhat smaller 
in magnitude. One also sees that the thresholds for quark quark scattering have 
shifted to slightly lower values, due to the temperature dependence of the quark 
masses themselves. 

Quark quark scattering processes that involve at least one strange quark are 
shown in Fig. ^ as a function of ^/s for the temperatures T = 215 MeV and T = 
250 MeV. The peak like structure in the us us scattering cross section is a cutoff 
artifact. The cross sections are otherwise featureless. Increasing the temperature 
from T = 215 MeV to T = 250 MeV reduces the us us cross section slightly, 
while raising the ss ss by a small amount. 

Differential cross sections may also be calculated for each process. As an example, 
we show the differential cross section da/d cos 6 for the process ud —>■ ud as a function 
of cos 6 at ^/s = 1 GeV and temperature T = 250 MeV in Fig. ^. For comparison, 
we also include this quantity calculated to lowest order in perturbative QCD and 
using finite values of the quark masses as are indicated by the NJL model at this 
temperature. An explicit formulation of the perturbative QCD cross section with 
finite masses is given in Appendix ^ Note that, since these cross sections are 
infrared divergent at this level, we have introduced an effective gluon mass WE 



ml = 2Tia, [l + ^]T^ ^ (600 MeV)^ , (2.24) 
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as a regulator, with A'^^ = 3 being the number of flavors here and the QCD coupling 
strength taken ad hoc to be = 0.6. The perturbative QCD cross section, 
which has essentially the same form as the M0ller scattering cross section for e~^e~ 
scattering, displays a strong forward peak. In the limit rUg 0, this peak becomes 
a pole and gives rise to the Coulomb singularity well known from QED. The NJL 
cross section, on the other hand, also displays a preference for the forward direction. 
However, this maximum is not as pronounced. 

C. Quark Antiquark Scattering 

We now turn to a discussion of processes of the form qq' qq'. For these, one 
has seven independent processes out of a total possible number of fifteen, taking 
again isospin and charge conjugation symmetry into account. These processes are 
listed in Table 0. The number of independent processes could be further reduced 
by taking into account crossing symmetry, e.g. by regarding us us and uu ss 
as dependent processes. In this paper, however, we regard processes, which are 
related by crossing symmetry as independent, since crossing symmetry does not 
lead to numerically equal cross sections. The relevant Feynman diagrams are shown 
in Fig. ^. Here only s-channel and t-channel diagrams occur. The species of the 
exchanged mesons are listed in Table y for the independent processes. 

As before, we will derive general expressions for the cross sections, illustrating 
these with the example us — > us. 

1. Analytical Calculations 

The Feynman diagrams for qq scattering are shown in Fig. ^. Analogously to 
Eqs. (|2.15| ), ( p.l6| ), the transition amplitude is given by 

- iMs = 6c^^c2^c3,c4v{p2)Tu{pi) iV^{pi+p2) u{p3.)Tv{p4) 

+ 6c^^c2^cz,cJ){p2){il^T)u{pi) [iV'^{pi+P2)\ u{P3)i^l5T)v{p4) (2.25) 

- iMt = Sc^^c3Sc2,cMp3)Tu{pi) [i'Cf (pi -ps)] v{p2)Tv{pi) 

+ Sc^,c35c2,cMP3)iil5T)u{pi) iV^{pi-p^) v{p2)in5T)v{p4) . (2.26) 
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The squares of these amphtudes can again be found in Appendix 0. 

We illustrate the further calculation by choosing the specific process us us. 
As can be seen from Table |T[ the s channel of this process proceeds via the exchange 
of a kaon in the pseudoscalar and a ax in the scalar part. This gives 

TPfT = 2Vk (2.27) 
TPfT = 2V„^ . (2.28) 

The t channel proceeds via the exchange of an r] or rj' in the pseudoscalar, a a or a' 
in the scalar part. This means that we have 



and an analogous expression for the scalar part. The differential cross section is now 

J2\Ms-Mtf (2.30) 



da 1 1 ■s-^ , . , . , |2 



dt 167r[s — {rriu + ms)'^][s — {rriu — rrisY] 
and the total cross section is calculated as 

a = Jdt^[l- iPE,)] [1 - fp i(3E,)] , (2.31) 

where a Fermi blocking factor for the final states has been introduced. In Eq. ( p.31| ), 

Ef = pf + mj, where i = 3,4. 



2. Numerical Results 



In Fig. 1^, we show a comparison of the scattering cross sections calculated in 
SUf{3) with the corresponding calculation in SUf{2), for a specific process ud ud, 
at the temperature T = 215 MeV. In contrast to the comparison shown in Fig. ^for 
quark quark scattering, the difference, on this scale, is not so great. This is due to the 
fact that the qq scatterings are resonance dominated: The temperature chosen lies 
only slightly higher than the pionic and kaonic Mott temperatures Tm-k = 212 MeV 
and Tmk = 210 MeV, so that both pions and kaons appear as very sharp resonances 
in the reactions containing them, see Table We therefore find that large cross 
sections occur at threshold for the process shown, in this case due to the pion 
resonance. We next show all qq cross sections that contain no strange quarks or 
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antiquarks in the incoming channel as a function of y/s at the temperature T = 
215 MeV in Fig. p. These cross sections all display pronounced structures, especially 
in comparison with the quark quark scatterings of similar flavor (cf. Fig. Q), as they 
admit resonances in the s-channel. Cross sections with uu in the initial state are 
in addition enhanced by the presence of a mesons. By contrast, the uu ss 
reaction contains no resonant structure and remains of the order of 1 mb over its 
range in ^/s. At the Mott temperature itself, intermediate states in the s-channel 
give rise to infinite cross sections at threshold. This feature, which also appears 
in other processes like vrvr tttc |T^, 717 iT'-f or qq 77 [^, is akin to 



the phenomenon of critical opalescence and has been discussed in greater detail 
in Ref. 0. For temperatures above the Mott transition, such as the one shown 
in this figure, the cross sections are large but finite at threshold. Increasing the 
temperature further to T = 250 MeV, as has been done for the same set of reactions 
in Fig. ^, shows that the resonances have become broader and the cross sections in 
this regime smaller. Nevertheless, this melting of the pion and sigma resonances still 
leads to cross sections that are highly enhanced over a range of ^/s between 0.2 and 
0.8 GeV in comparison with other scattering processes, where cross sections of only 
0{1 mb) are found. This remnant of the critical scattering phenomenon is thus still 
strongly visible even at high temperatures. Although it is not possible to observe 
such scatterings directly in experiment, it may lead to observable consequences when 
embedded in a consistently constructed chiral transport theory. 



In Fig. |T0|, we show the scattering cross sections for quarks and antiquarks for 
processes that contain at least one strange particle in the incoming channel at T = 
215 MeV. The threshold behavior is strongly dominated by the resonant structure 
for the case of the process us —>■ us. Here an exchanged kaon gives the dominant 
resonance, with a shoulder due to the ax- The ss uu singularity seen is a 
kinematical singularity, that arises because the reaction is exothermic. The ss ss 
displays no pronounced resonance - there is a small one due to the a' or rj'. At a 
somewhat higher temperature, T = 250 MeV, the resonant structure that forms a 
remnant of critical scattering for us — us, has become smaller and broader. Once 
again, it still peaks at 15 mb over a range of ^/s between 0.4 and 1 GeV, and may 
therefore also lead to an observable effect when embedded in a consistent transport 
theory. 
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Figure shows the differential cross section for the process us — > us, given as 
function of cos 6. Whereas the one gluon exchange perturbative QCD cross section, 
see Appendix Q once again resembles M0ller scattering for e^e~ and is strongly 
peaked in the forward direction, the NJL cross section is flat. This comes about 
due to the dominance of the s-channel exchange at this energy, which proceeds via 
spinless particles and thus shows no anisotropy. 

III. TRANSPORT PROPERTIES 

A. Averaged Transition Rates 

To calculate the transport coefficients, we need averaged transition rates. These 
are expressed as |^ 

in terms of the transition rate w{s,T), which in turn is defined as the product of 
cross section and relative velocity: 

w{s,T) = \vreiW{s,T) (3.2) 



\Vrel\ = ^ . (3.3) 

The quark density in Eq. (|3.1|) is defined as the integral of the Fermi distribution 
function over all momenta: 



To evaluate Eq. ( |3.1| ), we make the approximation that the total cross section is only 
a function of s, even if the incoming pair is moving with respect to the medium. In 
this case the averaged transition rate can be expressed as 

w{T) = f dsJ {piP2f - (mim2)V(s, T)P(s, T) , (3.5) 

J(mi+m2)^ 

where the weight function P{s,T) is given as 
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("OO 

x/ ciE2[2iVJ^(/3E2)]e(4|pinP2|'-(s-K + m2)-2EiE2)') . 
This weight function is illustrated in Fig. O for mi = 1712 = rriu, and the two 



values of the temperature T = 215 MeV and T = 250 MeV. One sees that this 
function weights the lower energies strongly and has an exponential tail at a given 
temperature. Increasing the temperature from T = 215 MeV to T = 250 MeV has 
the consequence that the peak value of the weight function is shifted and reduced 
somewhat and the exponential decay is weaker. Note that this approach differs from 
that taken in Ref. in that P{s,T) in Eq. (|3.6|) is not a normalized probability 



distribution. In Fig. 14, we show the averaged transition rates for the quark quark 
scattering processes uu —>■ uu and ud ud together with the quark antiquark scat- 
tering processes ud — > ud and us — > us. We see that the quark antiquark scattering 
averaged rates lie higher than the quark quark scattering ones. There is a sudden 
rise in the quark antiquark rates as one moves through the Mott temperatures. This 
is due to the inclusion of resonance channels at and above this point. Above the 
Mott temperatures, in the region of interest, both are decreasing with temperature. 

B. Thermal Relaxation Times 

Since the definition of the averaged transition rate strongly resembles the collision 
integral of a Boltzmann equation we can immediately identify the thermal 



relaxation time for each species as 

^7^ = Ep/'^//' ' (3-7) 
/' 

where the summation runs over all quark species and Wfg is the sum of the transi- 
tion rates of all processes with species / and g in the initial state. Our numerical 
calculation for = and is shown as a function of temperature in Fig. |15[ The 
temperature region below Tm-k is of purely academic interest, since physically there 
should be no free quarks: we show it here for completeness. The relaxation time is 
large for small temperatures, which is due to the low quark density. In the physically 
interesting region T > Tmtt, we obtain relaxation times of the order of 1-1.3 fm/c 
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for light quarks and 1.3-2 fm/c for strange quarks for temperatures greater than the 
pion Mott temperature. 



C. Mean Free Path 



We define the mean free path of a particle as [53 



A 



where the mean velocity of flavor / is given as 

2Nr f d^p p 



Our numerical calculation for A is shown in Fig. Once again, we also show the 
low temperature region, although, as already pointed out, this is unphysical; one 
can therefore regard the physical mean free path as being the value obtained in the 
high temperature regime. For high temperatures, T 3> mj, one finds Xf ~ Tf, since 
Vf is approximately equal to one for all flavors, as can be seen from Eq. ( |3.9| ). For 
lower temperatures, Vs < Vu < I, causing the curves in Fig. ITBIto lie closer together 



(3.8) 



(3.9) 



than the corresponding two curves of Fig. |I5|. From Fig. |I6| we obtain a mean free 
path of 0.9-1.4 fm for light quarks and 1.1-1.6 fm for strange quarks. 



D. Viscosity 

To a first approximation, the shear viscosity is proportional to the mean free 
path 

4 



Vf 
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PfPfX 



f 



(3.10) 



where P/ is the mean momentum for a quark of flavor /: 



Pf 



2iVe 
Pf 



d^p 
{2^ 



(3.11) 



A more detailed analysis p3| shows that in averaging the cross sections, large scat- 
tering angles should be preferred, which is achieved by replacing the total cross 
section implicitly contained in Eq. ( |3.10| ) with the weighted average 
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ldn^sm'e[l-Mm)][l-fFim)] . (3.12) 



Since in a first approximation, rj should be proportional to T^, we show r]u/T^ = 
rjfi/T^ and rjs/T^ in Fig. |l^. Since the mean free path for all flavors is approximately 
equal, the difference in the curves is mainly a density effect. Beyond the Mott 
transition, the result for rj/T^ lies in the range 0.83-1.0 for light quarks and 0.68- 
0.75 for strange quarks. The total viscosity, calculated for two flavors, amounts 
to 

3.2 < ri/T^ < 4.0 (3.13) 

and 

4.6 < ri/T^ < 5.5 (3.14) 

for three flavors. The value for two flavors lies under that given in Ref. 0. This 
is a direct manifestation of the difference in the cross sections shown in Fig. ^ for 
two and three flavors, due to the additional channels that are available in the three 
flavor calculation. A classical hydrodynamical estimate [0] places loose bounds on 
this number, 

2 < ri/T^ < 3tT , (3.15) 

in which r is the proper time for the expansion of the system. For t ^ 5 fm/c, 
our values for the total viscosity lie within these bounds, but are smaller than than 
those calculated in this hydrodynamical approach, which gives t]/T^ ^ 10 ||T2|. 



Alternatively, perturbative QCD calculations have been performed ||10| , pl]| for this 
quantity. The variational perturbative calculation of places a lower bound on rj 
as being 

/^•^ 1-16 

V/T'>~—, . 3.16 

aj log as 

For as = 0.6, this is 6.3, which is slightly higher than our result. Technically this 
comes about from the fact that the perturbative QCD calculation favors forward 
scattering, while the NJL model allows for a softer angular distribution. The per- 
turbative QCD result may however be substantially altered in the non-perturbative 
region, so we do not regard this as being a serious discrepancy. 
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IV. STRANGENESS ENHANCEMENT 



Given the elastic scattering cross sections, we are able to address the problem 
of strangeness enhancement in a quark meson plasma within our model. We choose 
an approach similar to that taken in Refs. p!3|j25|j26[] . However, instead of discussing 
gluon exchange, we compute the contribution of the exchange of mesonic resonances. 

The strangeness changing processes, which we have considered in this work, are 

uu ^ ss dd —>■ ss (4.1) 

and the respective back reactions 

ss uu ss —>■ dd . (4.2) 

In the following, we consider the light quarks to be fully equilibrated. We discuss 
two possible definitions of the chemical relaxation time. The first one is defined 
as the number of strange quarks present in chemical equilibrium, divided by the 
number of strange quarks produced per unit time: 

ri = ^^^ . (4.3) 

The factor 2 in this expression accounts for the number of light quark flavors, which 



give equal contributions. This quantity is comparable with that given in Ref. [|T3 
In this case one has, as a first approximation, 

p,(t)=p:'^tanh(t/ri) . (4.4) 

The second definition of the strange quark relaxation time is 

Again, the factor 2 counts the light quark flavors. This corresponds to making the 
rate equation ansatz 

-TT = —Pu , (4.6) 

dt T2 

neglecting the back reaction. The two definitions of Ti and T2 differ by the factor 
pf^{T) / pu{T) , which is smaller than one since ms{T) > mu{T), but approaches one 
for temperatures T ^ mg. 
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By contrast, in Ref. [0, the strangeness relaxation time was found within the 
framework of perturbative QCD to be dominated by the process gg — > ss. For this 
process, these authors obtained 

_ 1.61 (m,/T)exp(m,/r) 

a^T (1 + ||T/m, + . . .) ' ^^''^ 

where as is the strong couphng constant, which was assumed to be 0.6 for all tem- 



peratures. The relaxation time due to the processes (|4. 1| ) was calculated by the 
same authors to be roughly a factor of 10 larger. 

Our numerical calculation of ri, T2 and Tgg from Eqs. ( [4. 3D , (|4.5| ) and ( [4.7D are 



shown in Fig. |T8|. Of interest in this figure are the curves in the temperature range 
T > Tmtt- The solid curve shows ti, which lies in the range 11-17 fm/c, whereas T2, 
indicated by the dashed line, lies in the range 16-30 fm/c. The relaxation time due 
to gluon-gluon collisions, as given in Ref. |0 , is given by the dot-dashed line and is 
clearly seen to be considerably smaller than ti. However, according to Ref. [0, this 
calculation assumes a temperature independent value for the strange quark mass 
of rus = 150 MeV. The NJL model, on the other hand, predicts a strange quark 
mass, which, even at high temperatures, is at least twice as large as the current 
quark mass value of 150 MeV In the temperature range shown in the figure, 
150 MeV < T < 250 MeV, the strange quark mass falls from 510 MeV > > 
380 MeV. At the pion Mott temperature Tm^ = 212 MeV, we have = 415 MeV. 
An evaluation of Eq. ( ^4.7| ) with masses according to the NJL gap equation ( p.7| ) 
is indicated by the dotted curve in Fig. 113. Since the chemical relaxation time 
in Eq. (|4.7| ) depends exponentially on rris, an exact knowledge of this quantity is 
essential, as can be seen by comparing the dot-dashed and the dotted curve in 



Fig. |T8|. Taking the strange quark mass from the NJL gap equation gives a Tgg of 
the order of ti. Furthermore, Eq. ( [4.7| ) relies on the assumption, that light quarks 
and gluons are massless above the phase transition. Finite light quark and gluon 
masses could also shift the chemical relaxation times. 

Experimentally, the strangeness content of the observed mesons is parametrized 
in the ratio 

r (48) 
" 3 (it-) ' * ' 

where (Kg) and (vr") are the mean multiplicities of the observed mesons. At 

200 GeV/nucleon the experimental values are 
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(8.7 ±1.3)% N + N 

(4.9) 

(15.4 ±2.6)% S + S 



These values should be compared with the thermal value 

nh = ,J^f5^ (T^ = 27^^° (4.10) 
Pu{T)+pd{T) 

at T = 200 MeV. If we interpret the increase from r^"*"^ to r^^^ as being due to an 
"chemical equilibration" within the quark plasma alone, we can determine the time 
to available for equilibration from the equation 

tanh ( -) = ^^^^^ , (4.11) 

Vri/ rth 

and one obtains to ~ ti/4 ~ 4 fm/c. This value is clearly an upper limit, since 
strangeness can also be created within the hadronization process (found to contribute 
around 1% in our previous calculation [Q) and in the hadronic gas. This would mean 
that one has to assume a plasma lifetime of less than 4 fm/c. Note that calculations 
using the hydrodynamical Landau model [^, indicate that the total lifetime of the 
fireball (plasma ± hadronization ± hadron gas) is of the order of 5 fm/c for S + S 
collisions. This is not in contradiction with our estimate. 



V. GENERAL DISCUSSION 



In this section, we discuss general features relating to our calculation and its 
philosophy. We have evaluated the quark self energy in the Hartree approximation, 
which corresponds to the leading order term in a l/N^ approximation. The l/N^ 
expansion is necessary in order to have a well-defined method of dealing with a 
strong coupling theory. Concomitantly, the pion masses are calculated in the random 
phase approximation (RPA), which, taken together with the quark self energy, can 
be shown to form a consistent expansion in a chiral sense - i.e. the Goldstone 
mode is guaranteed in the chiral limit. The masses are not shown explicitly here 
- the calculation is standard and has been given in Refs. |^,|,^ in detail. In 
particular, to this level of approximation, the quark self energy is always real, even 
in a thermal medium, while the mesonic self energies are complex, depending on 
various parameters. 
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Given this circumstance, we have calculated quark-quark and quark-ant iquark 
scattering processes. In a simplistic fashion, by analogy with the Boltzmann equa- 
tion, we are able to extract transport coefficients, and in particular, lifetimes, which 
have been discussed in detail in Section |IT1[ These lifetimes imply of course that 
the u, d and s quarks carry widths that are inversely proportional to these. Thus 
one sees, that, while this is a consistent chiral and l/N^ expansion, it is not a self- 
consistent calculation of the widths. In order to do so, one must go beyond the 
lowest order l/N^ calculation. Several attempts to do this and also to include 
finite temperature [^] have difficulties either in including all possible terms, and/or 
maintaining chiral symmetry. A ffist calculation that attempts to go further 
and calculate the quark widths self-consistently ignores the consistent calculation 
of the meson sector, and thus violates chiral symmetry requirements at this level. 
We find that our widths ~ 1/r are of the order of 200 MeV and lie somewhat below 
the values quoted in Ref. PD . 



In this discussion, we wish to point out that there are contradicting expansions 
that one can examine. It has always been a fixed tenet of transport theory that 
truncations should be energy and symmetry conserving, and we have thus based 
our calculation on this. Clearly more work is required to incorporate the higher 
order terms in a self-consistent fashion that can also determine quark widths in 
medium simultaneously. As such, we can at best regard our transport theory results 
as estimates that may be compared with other similarly made estimates. 

Finally it is perhaps worth commenting that the correct inclusion of the quasi 



particles with widths is extremely important. In the calculation of Ref. |^], the 
authors have found that a consistent treatment of the widths leads to a correction 
of a factor of two for the collision rates. In addition, a consistent treatment may alter 
the transport coefficients substantially, when calculated field theoretically p2| , p3 



The imaginary part of the self energy can also be shown to be connected with 
memory effects . 



In concluding this section, it is evident that, within the NJL model, an expansive 
study that constructs a transport theory that accounts for all these problems, still 
requires substantial development. A ffist formal attempt at this has been made in 
Refs. |5|J5§. 
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VI. SUMMARY AND CONCLUSIONS 



We have examined the finite temperature behavior of elastic scattering processes 
of quarks with quarks and quarks with antiquarks within an effective three fia- 
vor model in which chiral symmetry is regarded as the most important underlying 
feature. In the quark quark elastic scattering, we have seen that the three flavor 
calculation differs quite substantially from a calculation that accounts for two flavors 
only, due to the influence of additional channels. The quark quark scattering cross 
sections are enhanced by a factor of 3-4, but are otherwise featureless. On the other 
hand the quark antiquark scattering cross sections display divergencies at the Mott 
temperature, which we can regard as an indication of critical scattering. At tem- 
peratures between Tmtt and T = 250 MeV, the remnants of the critical scattering, 
i.e. resonant structures that occur in the s-channel, are clearly visible and enhance 
the cross sections substantially in comparison with the background. Such a feature, 
it is hoped, may be visible in a final analysis that involves a consistent transport 
theory. The differential cross sections that we have calculated are found to be of the 
same order of magnitude as perturbative QCD calculations, in which a gluon mass 
of nig — 600 MeV and a strong coupling constant of ag — 0.6 have been assumed. 

We have also examined transport coefficients in the SUf{3) model. We find 
thermal relaxation times of 1-1.3 fm/c for light quarks and 1.3-2 fm/c for strange 
quarks. The mean free paths he in similar ranges, being 0.9-1.4 fm and 1.1-1.6 fm for 
light and strange quarks respectively. The ratios of the viscosity to the third power 
of the temperature, rj/T^ per quark fiavor, fall between 0.83-1.0 and 0.68-0.75 for 
light and strange quarks respectively. The total viscosity coefficient for two fiavors is 
found to be a factor roughly two smaller than previous SUf{2) calculations. This is 
attributable to the enhanced cross sections that enter into this SUf{3) calculation. 
Our value still falls within the bracketed range given by hydrodynamic estimates, 
while being slightly less than the perturbative QCD result. 

In studying strangeness enhancement due to the possible strangeness changing 
reactions that occur via quark antiquark annihilation, wc find that relaxation times 
are considerably larger than those derived from perturbative QCD with the naive 
assumption rris = 150 MeV. In the NJL model, nis = ms{T) is a function of tem- 
perature and its value lies significantly higher than the current quark mass value 
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rris ^ 150 MeV over the entire temperature range of interest. This feature follows 
from earlier temperature studies of the dynamically generated quark masses If 
one recalculates the perturbative QCD expression using the temperature dependent 
NJL quark masses, one arrives at chemical relaxation times that are of the same 
order as our calculation. In attempting to assess the relevance of elastic scatter- 
ing processes on strangeness production in heavy ion collisions, we note that this 
is a sensitive function of the plasma lifetime. A plasma lifetime of the order of 
4 fm/c would allow one to account for all of the observed strangeness enhancement 
in 5* + 5* collisions via the mechanism of elastic scattering. Since the plasma lifetime 
probably lies below this value, other processes such as hadronization and final state 
interactions must also play a role. 
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APPENDIX A: SQUARED TRANSITION AMPLITUDES 



1. Quark Quark Scattering 



The matrix elements in Eqs. ( p. 15 ) and ( 2.16 ) have the form 



iMt = Sc-,,c3Sc2,cMP3)u{Pi) il^hPi-Ps) u{pa)u{p2) 

+ Scuc-iSc2,c4u{Ps){il5)u{pi) iV[{pi-p3) u{p4){i-f5)u{p: 



(Al) 



and 



- iMu = Sc^^cJc2,C3u{Pi)uiPl) i'^uiPl-Pi) U{p^)u{p2) 

+ ^cr,ci^c2,c3u{pi){i-i^)u{pi) \iV^{pi -ps)! u{p:i){i-i^)u{p2) 



(A2) 



Here we have dropped the T factors for simplicity, since they can easily be included 
by a rescaling of the propagators. After a short calculation, one obtains 



25 



c s.c 



1-^*1' = ^{l^^ff tr[(^3+m3)(^i+mi)]tr[(:^4+m4)(#2+m2)] (A3) 
tr[(^3+ m3)75(^i+ mi)75]tr[(^4+ 7714)75(^2+ "^2)75]} ■ 



The interference term between the scalar and pseudoscalar part vanishes. The final 
result can be expressed using the Mandelstam variable t: 



1 



4A^2 

c s,c 



''13''24 



^13^24 5 



where we have abbreviated tf, = t— {rrii ±mjY. Analogously, one obtains 



4Ar2 

' c s,c 



"^14^23 + 



V. 



^14'"23 ' 



(A4) 



(A5) 



with Mj:- 



[mj im^)^. For the interference term, the spin summation results in 



4^EM-M: = 4^{PfPftr[(^3 + 



m3){pi+ mi)(p4+ m4){p2+ ^2) 



(A6) 



c 

- r'f'E'f*tr[(^3+ 7713) (|(i+ 7711)75(^4+ 7n4)(|*2+ 7772)75] 

- r'fDf*tr[(^3+ 7773)75(1*1+ 777i) (^4+ 7774)75(1(2+ 77I2)] 

+ A^^f*tr[(^3+ 7773)75(^1+7771)75(^4+ 7774)75(1*2+ "^2)75]} 

('^13^24 " ■5l2'^34 + "^li^tzj 



14 "23 



AN, 

— T)^T>^ (^13^24 ~ ^12^Zi + 

~ -^f (^^13^24 ~ '^12'^34 + "^14 "^23 

+ (^13^24 ^ '^12'^34 + "^14^23 



Here, we have abbreviated sfj = s— {rrii ± 777^)^. If the masses of all incoming and 
outgoing particles are equal, Eqs. (|A4|) -(|A^) can be largely simplified to yield 
1 



4Ar2 

c s,c 



EIMI 



Am 



2 _ 




^ {t - 4777^)2 + 




2 9 


(A7) 


2 _ 




^ {u - Arn^f + 


K 


2 2 


(A8) 



and 



^EM-M: = -^{Pfl)f[t« + 4777^(^ + 7/) - W] (A9) 

- w(t - 4777^) - PfPf t(« - Avr?) + PfPrtTz} , 

which is the result of Ref . H . 
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2. Quark Antiquark Scattering 



The matrix elements for qq scattering are 

- iMs = 4i,c2'^c3,c4^(P2)m(Pi) \i1^s (Pi + P2)] u{pz)v{pa) (AlO) 

+ ^c^,C2^C3,C4^v{p2){il^)u{pi) \iV^{pi+p2) \ U{p3){il5)v{p4) 



and 



- iMt = Sc^^c3^c2,cMPi)uiPl) [i'^tiPl -P3)\ V{P3)V{P2) 

+ Sci,c3Sc2,cMP4){il5)u{pi) \iVj^{pi -ps)] v{p3){i-f5)v{p 



(All) 



The square of these matrix elements can be immediately written down by using the 
fact that the crossing transformation 



-P4 



(A12) 



Pi -^Pl P2- 
P3 -P2 Pi P3 

transforms the t (u) channel of qq scattering to the s (t) channel for qq scattering. 
Thus one obtains 



4iV2 



4iV2 

' c s,c 



■^12-^34: + 



''13''24 ^ 



12 ■'34 



^13^24 1 



and 



4iV2 

c s,c 



E MsM: = — - [Pf (s+ S3+ - + t+3t+ 



(A13) 
(A14) 

(A15) 



■^12'^34 



'"l4'"24 + ^13^24 



"^f ('^12'^34 "^14^24 + '^13^24 



13''24 



Again, these expressions can be greatly simplified in the case of equal quark masses. 



APPENDIX B: ELASTIC SCATTERING IN PERTURBATIVE QCD 

In perturbative QCD, the elastic scattering of quarks and antiquarks proceeds in 
lowest order via one gluon exchange. A previous calculation for the case of massless 



quarks and gluons has been given in Ref . . Here we extend this result to the case 



of three flavors and introduce finite masses for quarks and gluons. 
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1. Quark Quark Scattering 



The elastic scattering of two quarks of flavor / and /' proceeds via a t-channel 
exchange, if / 7^ /', and via t- and ti-channel exchanges if / = /'. The differential 
cross section is 



da 



1 



1 



dt 167r[s — {rrif — mj/)^] [s — {rrif + '^/O^] 4Ai"^ 

For the t-channel exchange, one has, using the Feynman gauge 

1 



Y.\Mt-5fj>Mu\' . (Bl) 



iMt = u{p-i) ( -igs-i^-\% ] u{pi)- 



~-^^Hpd{-^9sYlxl)u{p,) , (B2) 



where we have introduced a finite gluon mass nig in order to avoid Coulomb singu- 
larities. For the -u-channel, one obtains 



u — mt 



(B3) 



Using standard techniques, the expression 



1 



c s,c 



2(s -mj- m),f + + 2st 



9 [ (t - m2)2 

2(s-2m^)2 + M2 + 2sM 



[u — m 



2^2 



2 (s - 4m|)2 - 4mj 



(t — ■m'i){u — m?) 



can be derived, where as = gl/^i-n) is the QCD fine structure constant. 



(B4) 



2. Quark Antiquark Scattering 

In the case of different incoming flavors, elastic scattering proceeds via a t- 
channel exchange. In this case, we have 

da 1 1 ■s-^ , . , i2 



dt 167r[s — {rrif — m^/)^] [s — {rrif + rrif/y] 4N^ 



EiMr , (B5) 



and the square of the scattering amplitude can be obtained from the t-channel 
amplitude contained in Eq. ((Bl]) by swapping s and u: 
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c s,c 



2^,2 



9" 



2(u 



m. 



m 



2 \2 



+ t^ + 2ut 



(B6) 



In the case of equal incoming flavors, we have to distinguish two possibilities: (i) the 
incoming and outgoing pair have equal flavor, and (ii) the incoming and outgoing 
pair have different flavor. Whereas the second process proceeds via an s-channel 
only, the first one also has a t-channel available. We now switch our notation to 
denote the flavor of the incoming pair by / and the flavor of the outgoing pair by 
/'. The differential cross section then takes the form 

da _ 1 1 

m ~ 167rs(s-4m/)2 4Ar2 



(B7) 



The squared transition amplitude can be obtained from Eq. (|B4| ) by making the 
substitution s ^ u ^ t ^ s, i.e. one has 



I 



4Ar2 

c s.c 



E \Ms - Sfj'Mt 



647r2a2 



In the massless limit, Eqs. 
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2(u 



m 



f 



m 



2 \2 



+ + 2su 



m 



2\2 



fB8) 



+ 5 



2(u 



2m)f + e 



9' 

2tu 



4mj 



3^'^' (s-mp(t-m2) 



(pq ) and (pq ) reduce to the result of Ref. p7 \. 



29 



REFERENCES 

[1] U. Heinz, Phys. Rev. Lett. 51 (1983) 351; Ann. Phys. (NY) 161 (1985) 48. 

[2] H. Th. Elze, M. Gyulassy and D. Vasak, Nucl. Phys. B 276 (1986) 706; Phys. 
Lett. B 177 (1986) 402. 

[3] U. Vogl and W. Weise, Prog. Part. NucL Phys. 27 (1991) 195. 

[4] S. P. Klevansky, Rev. Mod. Phys. 64 (1992) 649. 

[5] T. Hatsuda and T. Kunihiro, Phys. Rep. 247 (1994) 221. 

[6] J. Hiifner, S. P. Klevansky and P. Rehberg, University of Heidelberg Report 
No. HD-TVP-96-03 (1996), to appear in Physics Reports. 

[7] E. Quack and S. P. Klevansky, Phys. Rev C 49 (1994) 3283; J. Miiller and 
S. P. Klevansky, Phys. Rev. C 50 (1994) 410; V. Dmitrasinovic, H. J. Schulze, 
R. Tegen and R. H. Lemmer, Ann. Phys. (NY) 238 (1995) 332. 

[8] P. Rehberg, S. P. Klevansky and J. Hiifner, Phys. Rev. C 53 (1996) 410. 

[9] P. Zhuang, J. Hiifner, S. P. Klevansky and L. Neise, Phys. Rev. D 51 (1995) 
3728. 

[10] G. Baym, H. Monien, C. J. Pethick and D. G. Ravenhall, Phys. Rev. Lett. 64 
(1990) 1867. 

[11] A. Hosoya and K. Kajantie, Nucl. Phys. B 250 (1985) 666. 

[12] P. Danielewicz and M. Gyulassy, Phys. Rev. D 31 (1985) 53. 

[13] P. Koch, B. Miiller and J. Rafelski, Phys. Rep. 142 (1986) 167. 

[14] H. W. Barz, B. L. Friman, J. Knoll and H. Schulz, Nucl. Phys. A 519 (1990) 
831. 

[15] J. Bolz, U. Ornik and R. M. Weiner, Phys. Rev. C 46 (1992) 2047. 
[16] E. V. Shuryak, Nucl. Phys. A 566 (1994) 559c. 

[17] P. Rehberg and S. P. Klevansky, University of Heidelberg Report No. HD-TVP- 
95-13 (1995), |liep-ph/951022li 

30 



[18] E. V. Shuryak, Phys. Rep. 61 (1980) 71. 

[19] E. Quack, P. Zhuang, Y. Kalinovsky, S. P. Klevansky and J. Hiifner, Phys. Lett. 
B 348 (1995) 1. 

[20] A. E. Dorokhov, M. K. Volkov, J. Hiifner, S. P. Klevansky and P. Rehberg, 
University of Heidelberg Report No. HD-TVP-95-12 (1995). 

[21] P. Rehberg, Yu. L. Kahnovsky and D. Blaschke, University of Rostock Report 
No. MPG-VT-UR 81/96, in preparation. 

[22] S. R. de Groot, W. A. van Leeuwen, Ch. G. van Weert, Relativistic Kinetic 
Theory (North Holland, Amsterdam, 1980). 

[23] F. Reif, Fundamentals of Statistical and Thermal Physics (McGraw-Hill, New 
York, 1965) 

[24] M. H. Thoma, Phys. Lett. B 269 (1991) 144. 

[25] J. Rafelski and B. Miiller, Phys. Rev. Lett. 48 (1982) 1066. 

[26] T. Matsui, B. Svetitsky and L. D. McLerran, Phys. Rev. D 34 (1986) 783. 

[27] T. Alber et al., Z. Phys. C 64 (1994) 195. 

[28] P. Zhuang, J. Hiifner and S. P. Klevansky, Nucl. Phys. A 576 (1994) 525. 
[29] P. Zhuang, Phys. Rev. C 51 (1995) 2256. 

[30] E. Quack and P. A. Henning, Phys. Rev. Lett. 75 (1995) 2811. 

[31] L. Neise, V. Bunakov and J. Hiifner, University of Heidelberg Report No. HD- 
TVP-94-3 (1994), unpublished. 

[32] A. Hosoya, M. Sakagami and M. Takao, Ann. Phys. (NY) 154 (1984) 229. 

[33] P. Danielewicz, Ann. Phys. (NY) 152 (1984) 239. 

[34] P. A. Henning, Phys. Rep. 253 (1995) 235. 

[35] W. Zhang and L. Wilets, Phys. Rev. C 45 (1992) 1900. 

[36] W. Florkowski, J. Hiifner, S. P. Klevansky and L. Neise, Ann. Phys. (NY) 245 
(1996) 445. 

31 



[37] R. Cutler and D. Sievers, Phys. Rev. D 17 (1978) 196. 



32 



TABLES 



TABLE I. Independent processes for qq scattering. 



Process 


Exchanged Mesons in t Channel 


Exchanged Mesons in u Channel 


uu — ^ uu 


TT, rj, r]', (j^, (T, a' 


TT, 77, 77', (J^, (T, a' 


ss — ^ ss 


V, V', cr, 


77, 77', cr, cr' 


ud — ^ ud 


TT, ?7, 77', (Ttt, (T, a' 


TT, (Ttt 


us — > us 


V, V', 


K, gk 


TABLE IL Independent processes for qq scattering. 


Process 


Exchanged Mesons in s Channel 


Exchanged Mesons in t Channel 


ud — ud 


TT, (Ttt 


TT, 77, 77', (Ttt, cr, cr' 


us — ^ us 


if, (TK 


77, 77', (7, cr' 


uu — ^ uu 


TT, ??, 77', cr^, cr, cr' 


TT, 77, 77', cr^, cr, cr' 


uu — ^ dd 


TT, 77, 77', cr^, cr, a' 


TT, Cr^r 


uu — ^ ss 


r), r)', a, a' 




ss — uu 


77, 77', a, a' 


K, ax 


ss — ^ ss 


f], J]', a, a' 
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FIGURES 



FIG. 1. Feynman diagrams for elastic qq scattering to lowest order in l/Nc- 

FIG. 2. Comparision of the cross section for uu uu for three flavors (solid line) with 
the SUf{2) result of Ref. g (dashed hne), at T = 215 MeV. 

FIG. 3. Total cross section for elastic qq scattering involving only light quarks, as a 
function of ^s, at T = 215MeV and T = 250MeV. 

FIG. 4. Total cross section for elastic qq scattering involving at least one strange quark, 
as a function of ^/s^ at T = 215MeV (solid lines) and T = 250MeV (dashed lines). 

FIG. 5. Differential cross section for ud ^ ud sXT = 250MeV and y/s = 1 GeV within 
the NJL model (solid line) and perturbative QCD (dashed line). 

FIG. 6. Feynman diagrams for elastic qq scattering to lowest order in l/Nc- 

FIG. 7. Gomparison of the cross section ud ud calculated as a function of y/s and 
at T = 215 MeV for both SUf{3) (solid line) and SUf{2) (dashed line). 

FIG. 8. Total cross section for qq scattering with only light quarks in the initial state, 
shown as a function of ^/s, at T = 215MeV. 

FIG. 9. Total cross section for qq scattering with only light quarks in the initial state, 
shown function of at T = 250MeV. 

FIG. 10. Total cross section for qq scattering with at least one strange quark in the 
initial state, shown as a function of y/s, at T = 215MeV. 

FIG. 11. Total cross section for qq scattering with at least one strange quark in the 
initial state, shown as a function of y^, at T = 250MeV. 
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FIG. 12. Differential cross section for us ^ us at T = 250MeV and ^/s = 1 GeV 
within the NJL model (solid line) and perturbative QCD (dashed line). 

FIG. 13. Weight function P{s, T) as a function of s for mi = m2 = ruu, T = 215 MeV 
(sohd line) and T = 250 MeV (dashed line). 

FIG. 14. Energy averaged transition rates, shown as a function of the temperature. 
The pion Mott temperature is indicated by the dashed vertical line. 

FIG. 15. Thermal relaxation times for light (solid line) and strange (dashed line) 
quarks, shown as a function of the temperature. The pion Mott temperature is indicated 
by the dashed vertical line. 

FIG. 16. Mean free path for light (solid line) and strange (dashed line) quarks, shown 
as a function of the temperature. The pion Mott temperature is indicated by the dashed 
vertical line. 

FIG. 17. Viscosity coefficient per flavor for light (solid line) and strange (dashed line) 
quarks, shown as a function of the temperature. The vertical dashed line indicates the 
pion Mott temperature. 

FIG. 18. Chemical relaxation times for strange quarks, shown as a function of the tem- 
perature. The vertical dashed line indicates the pion Mott temperature. The relaxation 
time calculated according to Eq. ( |4.3| ) is given by the solid line, and according to Eq. ( [4.5D 
by the dashed line. The dot-dashed line gives the relaxation time calculated according to 
Eq. ( |4.7| ) with = 150 MeV, while the dotted line gives the relaxation time calculated 
according to Eq. ( |4.7] ), with strange quark mass from Eq. ( |2.7D . 
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